Transverse flow of nuclear matter in collisions of heavy nuclei at intermediate energies 
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The Quantum Molecular Dynamics Model (IQMD) is used to investigate the origin of the collective transverse 
velocity observed in heavy ion experiments. We find that there are three contributions to this effect: initial-final 
state correlations, potential interactions and collisions. For a given nuclear equation of state (eos) the increase 
of the transverse velocity with increasing beam energy is caused by the potential part. For a given beam energy 
the collective transverse velocity is independent of the nuclear eos but the relative contributions of potential and 
collisions differ. In view of the importance of the potential interactions between the nucleons it is not evident that 
the similarity of the radial velocities measured for fragments at beam energies below 1 AGeV and that for mesons 
at beam energies above 2 AGeV is more than accidental. 
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A major goal of heavy ion collision studies is to extend 
our knowledge of the properties of dense and hot nuclear 
matter. In the past these studies were concentrated on 
the production of secondary particles like pions or kaons 
and on anisotropies in the momentum distribution like 
the directed in-plane flow (bounce-off) or the out of plane 
flow (squeeze-out) Jl],|| of nuclear matter. For recent 
reviews of the different collective flows we refer to [^Jfl] 

Recently also the isotropic radial flow has gained large 
interest because at AGS and at CERN energies the ex- 
isting data at mid-rapidity of most of the mesons and 
baryons can be well described by the assumption that 
the system is thermalized and expands with a collective 
velocity |^|| . Later a similar observation has been made 
at energies between 150A MeV and 1A GeV (see eg. |?j 
- fig] ). There the fragments, observed at mid-rapidity, 
show an energy which increases linearly with the frag- 
ment mass. At these low energies most of the pions come 
from the decay of nuclear resonances and therefore can- 
not be used to determine a transverse velocity and kaons 
are too rare to be of use for a such an analysis. 

The radial flow velocity fit determined with help of 
the fragment kinetic energies up to 1A GeV smoothly 
joints that determined with help of mesons and baryons 
starting from 11 A GeV. This raised the question whether 
both transverse collective flows have a common origin 
and whether their functional dependence on the beam 
energy allows for conclusions about a possible transition 
between hadronic matter and a quark gluon plasma. 

A common origin of the transverse flow of fragments 
and mesons is not at all evident. Nucleons which are 
bound asymptotically in fragments do usually not suf- 
fer collisions with a large momentum transfer, otherwise 



they are not bound anymore. On the contrary, the cre- 
ation of mesons requires a large momentum transfer. 

At beam energies below lAGeV the radial flow seems 
to be a good candidate for measuring the nuclear equa- 
tion of state. Increasing compression leads to an increas- 
ing pressure which is released isotropically. The directed 
transverse flow in the reaction plane, which in principle 
also depends on the pressure, is known to be sensitive 
on the potential gradient. The squeeze-out of high en- 
ergy particles perpendicular to the reaction plane allows 
for a view inside the reaction zone. However, this effect 
is connected to the geometry and the time scales of the 
system. 

At CERN and AGS energies a detailed theoretical in- 
vestigation of the origin of the transverse velocity is dif- 
ficult jll| because many details of the expansion of the 
hadronic matter like the cross sections between baryonic 
and mesonic resonances are unknown. Therefore calcu- 
lations rely on a multitude of assumptions. 

At energies around 1 A GeV the situation is much bet- 
ter. Here we have models which simulate the complete 
reaction and we find agreement between the results of 
these models |IJ|^,[l2|-|l4| and experiment for many ob- 
servables. It is the aim of this article to take advan- 
tage of this situation and to investigate the origin of the 
transverse expansion of fragments in detail. Details of the 
Quantum molecular dynamics models, which are used for 
our study, can be found in the refs. |^,[l2| where also the 
different equations of state, mentioned in this article, are 
explained. We used for this work the IQMD version of 
the model with a variance of the Gaussian wave function 
in coordinate space of L — 4.33fm 2 and (if not otherwise 
stated) with static potentials because the compressional 
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energy is difficult to determine with momentum depen- 
dent interactions. A small value of L yields a higher 
number of fragments than a large one jl2],[l5| and repro- 
duces better the experimental multiplicity. This analysis 
is strongly based on the behavior of fragments. There- 
fore we have chosen this small L value. The choice of 
this value of L yield only a slight enhancement of the 
asymptotic value of (3 t for the 1 GeV case, where most of 
the analysis is done. For lower energies there is a visible 
reduction of the asymptotic (i t |1| . 

In the QMD simulations there are three different pro- 
cesses which contribute to the transverse velocity: 

• The nucleons which finally form a fragment A may 
have already initially a finite transverse velocity 



< Vi (A) >= ( - 




shows a linear velocity profile for central collisions at in- 
termediate energies. This linear profile is used as in- 
put to expansion models which assume the expansion 
to be due to a thermal and a collective part. Assum- 
ing E t (A) = l/2mAffi + 3/2T we disentangle these two 
parts by fitting the transverse velocities at mid-rapidity 
(y CM /Vpm\. ^ 0- 3 ) using two parameters fi t and T: 
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Am V Am 
where f3 t is the collective transverse velocity, T is the 
temperature, m is the mass of the nucleon and A the mass 
number of the fragments. v t (A) is the mean transverse 
velocity of fragments with the mass number A. We made 
two fits in order to be comparable with the two sets of 
experiments : one which includes all masses up to 
ten and a second where only masses smaller than five are 
included. 



where v\ equals vf (t — 0) and the sum runs over all 
nucleons which are finally entrained in the fragment 
of size A. <> means averaging over all fragments of 
size A observed at mid-rapidity (y CM / Z/p^j < 0.3). 
This velocity may depend on the size of the frag- 
ment. This type of initial -final state correlations 
have been investigated by Gossiaux et al. ]l7j . 

• The potential between the nucleons may contribute 
to the transverse velocity 
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To investigate its influence we sum up all contribu- 
tions from the nuclear force to the transverse ve- 
locity. 

• The collisions between the baryons may increase 
the transverse velocity as well 
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To calculate this contribution we add the trans- 
verse velocity change in all individual collisions 
(vj,V2) ~ > (v'jv'2) of the nucleons j of a fragment 
of the size A. 

Before we start with a detailed investigation of the 
origin of the transverse flow (3 t we compare our calcu- 
lations with the experimental results to make sure that 
the gross features of the transverse flow are well repro- 
duced in our simulations. As it was shown in |l5| ] in 
IQMD the transverse expansion of the nuclear matter 



Au+Au, very central collisions 



A IQMD b=lfm, fit for all A 

O IQMDb=lfm, fit for A <4 




• EOS: A <4, PRL 75 (95) 2662 

♦ FOPI: A <4, NPA 586 (95) 755 
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FIG. 1. Transverse velocity /3t obtained in IQMD simula- 
tions using a hard momentum dependent interaction as com- 
pared to the different experiments. 

Fig. 1 shows the results of the simulations for an im- 
pact parameter of b= 1 fm using a hard equation of state 
with momentum dependent interactions. These calcula- 
tions are compared with the available experimental data. 
The fit value of the different experiments have been deter- 
mined differently and in addition the impact parameter 
cut is different. The difference between the data of fij 
and of [|| is not understood yet. We observe that the 
simulations follow very closely the experimental data. 

We see as well in fig.l that even at the lowest beam 
energies a finite collective transverse velocity is observed. 
In this energy range Nebauer et al. discussed on cor- 
relations between the radial expansion and the directed 
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transverse flow. However, they studied semicentral col- 
lisions where the radial flow is small and the transverse 
flow is strong. In this article we will only regard very 
central collisions where the transverse flow is vanishing 
and the radial flow is strongest and analyze the influence 
of potential and collisions. 

Fig. 2 shows the transverse velocities taken at mid- 
rapidity for fragments produced in the reaction Au(250 
AMeV)+Au at b=lfm. The circles denote the final veloc- 
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FIG. 2. Transverse velocities vt(A) for the system Au(250 
AMeV)+Au. 

ities of the fragments. The dotted curve shows the result 
of a fit of the form of eq. 1 and illustrates the method of 
determining the asymptotic radial velocity as described 
above. 

The triangles mark the initial total transverse veloc- 
ity of the fragment nucleons. Clusters of nucleons which 
have a velocity pointing away from the interaction zone 
have a higher chance to survive the reaction without be- 
ing destroyed than those clusters which have a momen- 
tum which points into the reaction zone. For the dia- 
monds we have added to the initial transverse velocity the 
transverse velocity caused by the potential. The influence 
of the potential on the transverse velocity is weak at that 
energy. The difference between the diamonds and the 
circles, the total transverse velocity, is the transverse ve- 
locity caused by collisions. We see a continuous decrease 
of v t {A) with the fragment mass up to some asymptotic 
value. Thus, at 250 AMeV the transverse velocity has es- 
sentially two components: an initial transverse velocity 
and that caused by collisions. The potential plays only a 
minor role. 

It has now to be analyzed whether the vanishing in- 
fluence of the potential is a genuine effect or only valid 
for a beam energy of 250A MeV. Figure 3 shows the 
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FIG. 3. Energy dependence of the collective transverse ve- 
locities f3t for the system Au+Au and its different compo- 
nents. 

asymptotic velocity and its decomposition into the dif- 
ferent contributions. We display the initial transverse 
velocity of those nucleons which form finally a fragment 
as well as the sum of the initial transverse velocity and 
that caused by the potential. The difference between this 
curve and the asymptotic value is the transverse veloc- 
ity caused by collisions. For energies larger than 500 A 
MeV we see an almost constant contribution from the 
collisions ( of about (3 = 0.1) and from the initial ve- 
locity distribution of the fragments (of about (3 — 0.04). 
The increase of f3 t is almost exclusively caused by the 
potential. Higher beam energies yield a higher compres- 
sion and therefore stronger forces. Only at energies well 
above 2 GeV (where one may question the validity of 
our model) we observe saturation. The low energy sec- 
tor is more complicated. We see that below 250A MeV 
the forces become attractive (the reaction develops to- 
wards a scenario expected for deep inelastic collisions) 
and therefore they lower the contribution from the ini- 
tial velocity distribution. Collisions become more and 
more Pauli blocked and do not contribute anymore. At 
high beam energies the nucleons which are finally en- 
trained in fragments have been (in central collisions) ini- 
tially close to the surface at the back end of the nuclei, 
because there the number of collisions they suffer is low- 
est [fj"7f . Also their initial momentum points away from 
the reaction zone. These initial-final state correlations 
change at lower energies. There, due to the increased 
Pauli blocking, projectile matter can traverse the target 
and vice versa and hence the nucleons which are finally 
found in fragments have initially a quite uniform distri- 
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FIG. 4. Time evolution of the density at the center of the 
reaction zone for equations of state with different compress- 
ibilities. 

bution over the nucleus. In Jl6| an initial final state cor- 
relation is found for fragments emitted under 9cm = 90°. 

The physical origin of the energy dependence of (3t can 
also be studied by comparing it with the radial flow ex- 
pected if the total compressional energy is converted into 
a radial flow. For this purpose we compare the energy 
difference AU = U max ~ Um where U max is the total po- 
tential energy at the moment of the highest compression 
whereas Ui n i is the initial potential energy of projectile 
and target. The difference corresponds to the compres- 
sional energy needed to reach the state maximum com- 
pression. 

AU is related to /3 t by 
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The values are marked by diamonds. We see that its 
excitation function agrees to that of the contribution of 
initial (3 t plus potentials (triangles), which is the total 
radial flow minus the collisional part. For higher ener- 
gies this curve is parallel to that of the total radial flow 
(circles) what verifies that the increase of the radial flow 
is caused by the increasing compression, whereas the col- 
lisional contribution stays rather constant. 

All calculations discussed up to now have been per- 
formed using a hard equation of state with momentum 
dependent interactions (mdi). We now compare differ- 
ent equations of state at an energy of 1A GeV. In the 
Skyrme ansatz the repulsion of the nuclear equation of 
state at high densities is connected to the nuclear com- 
pressibility. Therefore we characterized in fig. 4 and 5 the 
equation of state by their corresponding compressibility. 
Soft and hard eos correspond to k = 200 and n = 376 
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FIG. 5. Dependence of the collective transverse velocities 
j3t on the compressibility for the system Au(lA GeV)+Au 
and the contribution of the potential part. 

MeV, respectively. Fig. 4 shows the time evolution of 
the density in the center of the reaction zone of a central 
Au+Au collision using different equations of state with- 
out momentum dependent interaction. We see that due 
to less repulsion the highest density is reached with the 
softest equation of state. Also the pressure is lower for 
this eos. Consequently the system will stay for a longer 
time at maximum compression. A higher density causes 
a smaller mean free path and therefore an enhanced num- 
ber of collision. Both yield together a nearly linear de- 
pendence of the number of high energy collisions on the 
compressibility. Independent of the eos each high ener- 
getic nn collision contributes about the same amount to 
the radial velocity. The softer the eos the more collisions 
contribute to the radial velocity and the higher is the 
kinetic pressure caused by the collisions. On the other 
hand the repulsion of the potentials is weaker and causes 
a weaker pressure from the potential part. Therefore it 
has to be checked which contribution to the pressure is 
dominant and whether the dependence of the different 
contributions to the total pressure counterbalance. 

Figure 5 shows the asymptotic transverse velocities 
and the contribution of the potential part as a func- 
tion of the compressibility for four different equations 
of state of Skyrme type. The difference is caused by col- 
lisions and the initial state collective velocity. The latter 
contribution may be regarded as independent to the nu- 
clear eos. The values for a hard eos (x=376 MeV) resp. 
soft (x=200 MeV) without momentum dependent inter- 
actions are very comparable to those with mdi. It is a 
remarkable result that all equations of state yield about 
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the same fit despite the contributions from the potential 
part are rather different. As already stated the potential 
part and the collisional part show different dependences 
on the nuclear eos. Seemingly the differences counterbal- 
ance. 

In fig. 5 we show as well as open triangles the radial ve- 
locity one would obtain if the compressional energy per 
nuclcon is entirely converted into a radial velocity apply- 
ing eq.2.: We see that this line parallels that with the 
actual final transverse velocities caused by the potential. 
Thus the increasing compressional energy is responsible 
for the increase of the radial velocity. The different abso- 
lute values are due to the fact that the system is not in 
complete equilibrium and hence the compression is low- 
ered due to the not completely decelerated projectile and 
target matter. 

Using the IQMD model to investigate the collective 
transverse expansion velocity observed in heavy ion re- 
actions we find, as in the experiments, a continuous rise 
with the beam energy up to energies of 1600A MeV. 
Whereas the in-plane flow measures the potential gra- 
dient and the out of plane squeeze the mean free path 
it was hoped that the radial flow is directly related to 
the compressibility. This conjecture can only partially 
be substantiated. Indeed, the increase of the radial flow 
as a function of the bombarding energy is - for a given 
nuclear equation of state - due to the higher compression 
of the nuclear matter and not due to the rescattering 
of the fragment nuclcons. The relative fraction of the 
potential and of the collisional contribution at a given 
energy depends on the other hand on the compressibil- 
ity of the equation of state. Even more, the sum of the 
potential and collisional contribution is almost constant. 
Thus it is impossible to extract the compressibility from 
the observation of the radial flow. 

The collisional contribution to the radial flow should be 
rather similar for the mesons as well as for the baryons. 
At energies around 1 AGeV, however, the potential con- 
tributes in between 30% and 60% to the total radial ve- 
locity, depending on the equation of state. This contri- 
bution is absent for the mesons if created by string decay 
and not by baryonic interactions like NN — ► A7V. Be- 
cause no potentials are at hand to continue our investi- 
gations to higher energies the doubts remain whether the 
one should expect a smooth transition between the collec- 
tive radial velocity observed 1 AGeV for fragments and 
that measured for mesons at higher (AGS and CERN) 
beam energies. 
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